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The superconducting state of Ca-intercalated graphite CaC6 has been investigated by specific 
heat measurements. The characteristic anomaly at the superconducting transition (T c = 11.4 K) 
indicates clearly the bulk nature of the superconductivity. The temperature and magnetic field 
dependence of the electronic specific heat are consistent with a fully-gapped superconducting order 
parameter. The estimated electron-phonon coupling constant is A = 0.60 - 0.74 suggesting that the 
relatively high T c of CaCg can be explained within the weak-coupling BCS approach. 



PACS numbers: 74.70.-b, 74.25. Bt, 74.70.Ad 



The recent discovery of superconductivity in Ca- and 
Yb-intercalated graphite has refocussed considerable in- 
terest onto graphite intercalated compounds (GICs) 0, 
2]. The superconducting transition temperature for Ca- 
and Yb-intercalated graphite is T c ks 11.5 K and 6.5 K, 
respectively, significantly higher than the alkali-metal in- 
tercalated graphite phases studied in the 1980's 3]. Sim- 
ilar to MgB 2 where the hexagonal B sheets are interca- 
lated with Mg, in the GICs the Ca or Yb atoms are sand- 
wiched by the honeycomb graphene layers. The interca- 
lated metal ions act as donors and transfer charge into the 
host graphene layers, resulting in partially filled graphene 
7r-bands. In contrast to MgB2 with a strong electron- 
phonon coupling of the B cr-bands, the coupling strength 
of the graphite 7r-bands to in-plane phonon modes is ex- 
pected to be small. Also out-of-plane phonon modes will 
not couple to the 7r-band because of the antisymmetric 
character of the 7r-orbitals. Thus the superconducting 
mechanism for the GICs could be rather different than 
that of MgB 2 . 

The superconducting mechanism, now under debate, 
relies on the other partially filled bands: the so-called in- 
terlayer bands. From the detailed comparison between 
the superconducting and non-superconducting GICs, 
Csanyi et al. found that the filling of the interlayer state 
is essential for stabilizing the superconductivity in GICs 
■ Based on this result and the consideration about the 
layered structure of the GICs, they proposed an uncon- 
ventional electronic paring mechanism involving excitons 
or low-energy acoustic plasmons 0. The latter has 
been discussed as a possible mechanism for the high T c of 
the intercalated layered metal halide nitrides Jb| . Soon af- 
ter, however, Mazin || and Calandra et al. yj suggested 
that the interlayer band originates from the s-band of 
the intercalant, and an ordinary electron-phonon cou- 
pling mechanism involving the intercalant phonon and 
the out-of-plane C phonon modes would be sufficient to 
explain the relatively high T c in the GICs. In order to 
shed light on this controversy, further experimental stud- 



ies on the superconducting properties of the GICs are 
required. 

Specific heat (Cp) measurements are known to be a 
powerful tool to investigate the superconducting states. 
Since Cp probes the quasi-particle excitation across the 
superconducting gap, its temperature and magnetic field 
dependence reflect the nature of the superconducting 
state such as gap symmetry, presence of multi-gaps, and 
coupling strength between electrons and phonons (or 
other bosonic excitations). So far, there has been no Cp 
study on the superconducting state of the GICs, mostly 
because the low T c was experimentally difficult to reach, 
and the sample homogeneity was not sufficient for reliable 
studies. In this Letter, we report the temperature and 
magnetic field dependence of Cp for high-quality bulk 
CaC6 samples. In zero-field measurements we observe a 
sharp Cp anomaly at T c ~ 11.4 K. The detailed temper- 
ature and magnetic field dependence of Cp clearly show 
a fully gapped superconducting order parameter, and fit 
very well to the weak coupling BCS prediction. 

The Ca-intercalated graphite samples were synthesized 
by reacting highly oriented pyrolytic graphite (HOPG, 
Alfa Aesar, spread of the c-axis orientation 0.4(1)°) with 
a molten lithium-calcium alloy at 350° for several weeks 
(cf. Ref. 0)- Since the metal alloy and the interca- 
lated sample are air sensitive all handling was done in 
purified argon atmospheres. X-ray diffraction patterns 
(not shown) dominantly show reflections of CaCe (the 
distance between graphene layers, d = 4.50 A), and con- 
tained no reflections related to non-intercalated graphite. 
Weak additional reflections jlCJl possibly due to Li/Ca- 
intercalated graphite phases [ll| were observed. The con- 
tribution of the impurity phase is estimated to be less 
than 5% to the total heat capacity. The magnetic sus- 
ceptibilities were determined in a MPMS SQUID magne- 
tometer (Quantum Design). The superconducting tran- 
sition temperature of our samples is T c — 11.40(5) K 
see the inset of Fig. 1), consistent with previous reports 
The transition width determined as the temperature 
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difference between 10% and 90% diamagnetic shielding is 
only 0.1K, smaller than found in previous reports [2| . Af- 
ter demagnetization correction, the volume fraction was 
estimated to be « 100 %. The heat capacity for a sam- 
ple of « 5 mg was measured using a PPMS calorime- 
ter (Quantum Design) employing the relaxation method. 
To thermally anchor the crystals to the sapphire sample 
platform and to orient the samples with respect to the 
magnetic field (H || c), a minute amount of Apiezon N 
grease was used. The heat capacity of the platform and 
the grease, which amounted to ~ 40-60% over the whole 
temperature range, was determined in a separate run and 
subtracted from the measurements for the samples. 

Figure 1 shows the temperature dependence of Cp at 
H = and H = 10 kOe applied along the c-axis. A sharp 
anomaly at T c is resolved indicating the bulk nature of 
the superconductivity in our sample. The onset of 11.40 
K determined from the specific heat jump is consistent 
with the T c obtained from the susceptibility measure- 
ments. A small offset of Cp/T at H— (7„ s = 0.65 ± 
0.13 mJ/mol K 2 ) and a slight hump of C P /T at H = 10 
kOe were observed below 3 K (see the right inset of Fig. 
1), which probably originate from paramagnetic impuri- 
ties or a non-superconducting fraction in our sample. 7„ s 
is only ~ 10 % of the Sommerfeld coefficient 7^ extracted 
from the normal state Cp(T) as discussed below. If we 
assume that 7„ s corresponds to the non-superconducting 
part of our sample, the volume fraction of the supercon- 
ducting portion can be estimated as 1— ( r y ns /'J%) ~ 90% 



proving a good quality of our GIC sample. The upper 
critical field along the c-axis has been reported to 
be ~ 2.5 kOe |2|], thus the normal state Cp(T) can be 
obtained when the superconductivity is completely sup- 
pressed by a magnetic field H = 10 kOe. A clear devi- 
ation from Debye T 3 law seen in the normal state Cp 
can be attributed to a contribution from the low-lying 
optical phonon modes. Recent a6-initio calculations also 
show the presence of low frequency modes mainly involv- 
ing vibration of the intercalated Ca . The normal state 
Cp can be described by Cp — ~f* N T + Ci a ttice{T) where 
7^T is the electronic contribution, and Ci a ttice(T) = /3T 3 
+ 5T 5 is the lattice contribution. The solid line in Fig. 1 
is the best fit to the H = 10 kOe data below T < 12 K, 
yielding the parameters, 7^ = 6.66(1) mJ/mol K 2 , (3 = 
65.1(3) AiJ/mol K 4 , and 5 = 0.256(3) /iJ/mol K 6 . Con- 
sidering the offset of Cp/T(H = 0), j na , the Sommerfeld 
coefficient for the superconducting part of the sample is 
7jv ~ l* N -lns = 6.01 ± 0.14 mJ/mol K 2 . The corre- 
sponding Debye temperature 0_d(O) is 598(3) K, which 
is much higher than that of pure graphite (6_d(0) = 413 
K) and the alkali metal GICs (e.g. 235 K for KC 8 ) fTl, 
but comparable with 613(0) = 590 - 710 K for LiC 6 [lg. 

The specific heat difference ACp between the normal 
and superconducting state is shown in Fig. 2. The solid 
line is the theoretical fit based on the 'a-model' assum- 
ing an isotropic s-wave BCS gap A(T) scaled by the ad- 
justable parameter, a = A(0) /kpT c 14]. For the weak 
coupling limit a is 1.76. The detailed temperature de- 
pendence of ACp/T was fitted by three adjustable pa- 
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FIG. 1: (Color online) Temperature dependence of the spe- 
cific heat at H= and 10 kOe. The solid (black) line through 
the data points at H = 10 kOe (also in the right inset) is a 
fit described in the text. (Left inset) The temperature de- 
pendence of the susceptibility with H — 6 Oe under zero-field 
cooled(ZFC) and field-cooled(FC) modes. (Right inset) The 
low temperature behavior of Cp(T). The (red) solid line for 
C P (T) at H = is a linear fit for 2 K < T < 3 K. 
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FIG. 2: (Color online) Temperature dependence of ACp/T 
= C P (H= 0)/T - C P (H = 10 kOe)/T. The (red) solid line 
is the best fit according to the a-model assuming an isotropic 
s-wave BCS gap as described in the text. The electronic con- 
tribution of the specific heat C es is plotted on a logarithmic 
scale against T c /T in the inset. The (red) straight line in the 
inset is the exponential fit for 2.5 < T c /T < 5.5. 
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rameters: a, Sommerfeld coefficient 7, and T c . The data 
are very well reproduced by a = 1.776, 7 = 5.91 mJ/mol 
K 2 and T c — 11.30 K. The normalized specific heat jump, 
ACp/'yT,, is 1.432, which is close to the weak limit BCS 
value, 1.426. The normalized electronic specific heat, 
C es /jT c in the superconducting state is shown in the in- 
set. The solid line is an exponential fit to the data for 2.5 
< T c /T < 5.6 using the form C es /jT c oc exp(-0.82<xT c /T) 
with a =1.65(1). C es exponentially vanishes for T — > 
K, clearly manifesting the absence of gap nodes in the su- 
perconducting state. The a value is somewhat lower than 
that from the a-model fit in which the ratio A(0)/A;bT c 
is mostly determined by the shape of the Cp jump near 
T c . However the discrepancy is less than 10%, in contrast 
to MgB 2 where the discrepancy is more than 60% due to 
the presence of two gaps [l5j . 

Figure 3 shows the temperature dependence of Cp/T 
measured with different magnetic fields. Superconduc- 
tivity is gradually suppressed, and the quasiparticle con- 
tribution increases with magnetic field. T C (H) deter- 
mined by the mid-point over the Cp anomaly is plotted 
in the inset with the error bar corresponding to the tran- 
sition width. The T C (H) obtained from the susceptibility 
measurements (not shown) are in good agreement with 
those found from Cp(T). For comparison, we also plot 
the predicted curve based on the Werthamer-Helfand- 
Hohenberg (WHH) theory 0. H^(T) shows a linear 
dependence down to T/T c « 0.2 and a clear deviation 
to the WHH prediction, which deserves further investi- 



gation. The best fit for low fields results in (dH^/dT)^ 
= 219.1(8) Oe/K and the linear extrapolation of H^(T) 
to T = K yields ff^(0) = 2.48 kOe, which corresponds 
to an a6-plane coherence length, £at>(0) ps 360 A. From 
i?| 2 (0) w 7 kOe along the a6-plane 0, the c-axis coher- 
ence length, £ c (0) ~ 130 A (> d = 4.50 A) is obtained. 
This result indicates the 3-dimcnsional nature of super- 
conductivity in C&Cq 0- 

The magnetic field dependence of the Sommerfeld co- 
efficient j(H) provides information about the supercon- 
ducting gap symmetry. For a gapped superconductor, 
the quasiparticle excitations are expected to be confined 
in the vortex cores, and 7(H) is proportional to the den- 
sity of vortices resulting in 7(H) oc H. In contrast, when 
the superconducting gap is highly anisotropic or has a 
gap node, the delocalized quasiparticles near the gap 
minima cause a nonlinear magnetic field dependence such 
as 7(H) oc H 1 / 2 [l^. The normalized 7(H) for CaC 6 is 
displayed in Fig. 4. 7(H) for T — > K is determined from 
the linear fit for 2 K < T < 4 K (see the inset of Fig. 
4) after subtracting the offset of 7 ns (H = 0). As a con- 
sistency check, we also plot C P {T,H)/T (» C es {Tfl)/T 
+ y(T,H) + Ciattice/T) at the lowest temperature in 
the present measurements (T = 2 K) after subtraction of 
C P (H = 0)/T(ra C es {Tfi)/T + C latUce /T). Both curves 
are normalized with y(H — 10 kOe). 

As shown in Fig. 4, 7(H) increases linearly with H for 
low fields up to H w 0.3 H C 2- This behavior is strongly 




FIG. 3: (Color online) Magnetic field dependence of the heat 
capacity for CaC6 . The numbers next to the data correspond 
to the applied magnetic field (kOe) along the c-axis. The inset 
shows the H^(T) for H \\ c estimated from specific heat (■) 
and susceptibility (O)- The (blue) dashed line demonstrates 
the WHH prediction Hql . and the (red) solid line is a linear 
fit for the low magnetic field data (H < 1 kOe). 



FIG. 4: (Color online) Magnetic field dependence of 7(H) 
as a function of 7(H) for T -> K (■) and at T 

= 2 K (O) is normalized by 7^ = 7(10 kOe). The (red) 
solid and (blue) dashed lines correspond to 7(H) oc H and 
7(H) oc H 1 / 2 , respectively. The (black) dotted line indicates 
y(H)/7N = 1. (Inset) Low temperature Cp data at various 
magnetic fields. The (black) solid lines are fits to the data for 
2 K < T < 4 K. The numbers next to the data correspond to 
the magnetic fields (kOe). 
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adverse to j(H) oc H 1 / 2 expected for nodal supercon- 
ductivity and also in contrast to two-band superconduc- 
tivity , v lj|. Recent calculations found that even for an 
isotropic-gapped type-II superconductor, j(B) oc B be- 
havior persists only up to a certain crossover field, B* 
[l9|. For the isotropic superconducting gap, B* is ex- 
pected to be w 0.32B c2 which is reduced as the degree 
of the anisotropy for the superconducting gap increases. 
Since Cp{H) measurements were done in the field-cooled 
mode, we can approximately estimate the ratio between 
the crossover field and the upper critical field B*/B C 2 
w H*/H C 2 ~ 0.3, which is quite close to the isotropic 
limit. Thus the magnetic field dependence of 7(H) con- 
sistently supports the notion of a fully-gapped and almost 
isotropic superconducting order parameter. 

Finally, we estimate the electron-phonon coupling 
strength, A for CaCg. At first, A can be estimated 
from a comparison between 7jv and the density of 
states at the Fermi level Ep, N(0) using the equation 
7JV = (27r 2 fc|/3)iV(0)(l+A). If we take N(0) = 1.50 
states/eV-cell from recent calculations 0, we obtain A « 
0.66 - 0.74 using our 7jv = 6.01 ± 0.14 mJ/mol K 2 . Alter- 
natively, we can also estimate A based on the McMillan 
formula 1201, 



O, 



1.45 



exp 



-1.04(1 + A) 
A - (l + 0.62A)/i* 



(1) 



Here, /i* is the Coulomb pseudopotential. Taking T c — 
11.4 K and 0d = 598(3) K from the normal state Cp, we 
obtain A = 0.60 - 0.71 with \i* ranging from 0.10 to 0.15. 
All these results are in good agreement with each other, 
and imply that CaCg obviously belongs to the weak cou- 
pling regime. Recent calculations [9j also predicted A to 
be w 0.83, in accordance with the measured A. Consis- 
tently, the superconducting gap ratio, 2A(0)/kpTc = 3.3 
- 3.6 is also close to the BCS value 3.52. 

In conclusion, we have reported the first specific heat 
measurements for the superconducting and normal state 
of high-quality bulk CaC6 samples. The specific heat 
anomaly at T c is clearly resolved indicating the bulk 
nature of the superconductivity. Both the tempera- 
ture and magnetic field dependence of Cp strongly sug- 
gests a fully-gapped superconducting order parameter. 
The estimated electron-phonon coupling constant A is 
0.60 - 0.74, in a good agreement with recent calcula- 
tions. These results suggest that CaC6 is a fully-gapped, 
weakly-coupled, phonon-mediated superconductor with- 
out essential contributions from alternative paring mech- 
anisms. We cannot rule out the possibility of a super- 
conducting gap anisotropy or a multi-gap scenario in this 
system |2l|, however, if any, the difference between the 
largest and smallest gap s will be rather small. Very re- 
cently, Lamura et aZ.[22j reported the in-plain penetra- 
tion depth, A a h for the bulk CaC6 sample. A clear expo- 
nential temperature dependence of A a & indicates a fully- 
gapped superconductivity in CaCg, which is consistent 



with our Cp study. 
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